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ABSTRACT 

We describe numerical methods for constructing lump solutions in open 
string field theory. According to Sen, these lumps represent lower dimensional 
Dp-Branes and numerical evaluation of their energy can be compared with 
the expected value for the tension. We take particular care of all higher 
derivative terms inherent in Witten's version of open string field theory. The 
importance of these terms for off shell phenomena is argued in the text. 
Detailed numerical calculations done for the case of general p brane show 
very good agreement with Sen's conjectured value. This gives credence to 
the conjecture itself and establishes further the usefulness of Witten's version 
of SFT . 



1 Introduction 



Since its introduction, string field theory held out the promise for nonper- 
turbative studies of string theory. Recently Sen [l], || has argued that open 
bosonic (string) field theory describes the dynamics of DD system with the 
tachyon providing the instability inherent in such pair. Tachyon condensation 
also describes the decay of a single unstable D-brane. In addition, according 
to @ the kink and lump solutions of such field theory lead to lower dimen- 
sional branes. In recent work Sen and Zweibach || studied in detail tachyon 
condensation in 26D open string theory. This follows the earlier, pioneering 
work of Samuel and Kostelecky who were the first to consider the vacuum 
structure of string field theory [|], |[. They used a level truncation scheme 
to generate an approximation for the tachyon effective potential. Following 
this scheme, results were found that show great agreement in numerical 
values with the expected exact results. Similar results were recently obtained 
for superstring field theory || [|J. Impressive high level studies appeared in 

SB- 

It is equally important to give a construction of non-constant kink and 
lump-like solutions. One can expect that it is for these that the stringy 
effects present in string field theory might play the most important role. 
One of the characteristic features present in construction of Witten's version 
of the theory 03, 0, |TJ, IH is the appearance in the interaction of terms 
exponential in derivatives. These terms can be moved from the interaction 
to give a nontrivial kinetic term. They have frustrated early attempts for 
construction of nonperturbative soliton (or instanton) solutions of the theory. 

In the present work, we consider this problem in numerical terms. Con- 
centrating on the lump of open string theory, we develop methods for its 
numerical solution. In this we keep the nontrivial exponential terms char- 
acteristic of the string vertex interaction and simultaneously perform a level 
truncation. This, we argue is to provide a very good approximation to the 
exact result. 

The content of the paper is as follows. After a short description of open 
SFT, we discuss some features relevant to the present work. We explain 
(based on earlier observations ||14|| ) how and why the approach of keeping 
higher derivatives and simultaneous level truncation holds the promise for a 
good approximation. We then proceed to the numerical work. 
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While this work was in progress, there appeared the work of Ref . [|T5f 
which considers the problem in its field theory limit. 



2 Open String Field Theory 

We begin by describing some features of open string field theory which are 
of relevance to the investigation that follows. One has the cubic action: 

S = (A\Q\A) + 9 -(V 3 \\A)\A)\A) (1) 

with Q being the first quantized BRST operator. The geometric, three string 
interaction is realized in the Hilbert space by the vertex 

(V 3 \ = (0|exp{£a/iO4 l + c/ln 7 £d r 2 } (2) 

r=l 

Here the Neumann coefficient N rs are determined in terms of appropriate 



conformal mapping , their explicit values are determined in PJ| . 

One of the main properties of the three-string vertex is an explicit ap- 
pearance of higher derivative terms. They come in exponential form acting 
on each string field 

\A) ^exp(a'hi7<9 M «^)|A) (3) 

with the constant 

1 = — (4) 

The exponential terms are not relevant in studies of vacuum structure but 
they can have a nontrivial effect in any other nonzero momentum process. 
Concerning a systematic approximation or expansion scheme one notes the 
following. The masses of tachyon (and higher mass fields) are proportional 
to 1/a! . In general a' serves as a scale ,it can be scaled out in front of the 
SFT action. At a nonperturbative level, there is no small free parameter 
and no systematic expansion. Since one is not able (at present) to solve the 
theory in exact terms, one relies on seemingly and hoc approximations. Such 
is the process of level truncation. In order to understand more clearly the 
procedure involved and the relative relevance of particular terms, let us recall 
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the original argument given for the level truncation in an unpublished work 



of ref.[14]. Considering an approximate calculation of a nonzero momentum 



amplitude, for example for four tachyons one starts from: 

A s = f 1 dxx- s / 2 - 2 (V u (3')\b x R \V 12 (3')), (5) 
J o 

representing the s-chanel Feynman diagram. 

Using the fact that well as similar results for b's and 

c's, we have 

A 3 = ^ dxx^ 2 - 2 (V M (3)b \V 12 (3)). (6) 
jo 

The dot on V indicates that a_ n , b_ n , c_ n have been replaced by x n 
a- n ,x n b- n ,x n C- n . Expanding in levels corresponds to expanding in pow- 
ers of x = e~ r . With the use the appropriate Neumann coefficients after a 
straightforward algebra one has 



A s = £ dxx- s ' 2 - 2 (l - l -^x 2 + ■ ■ e E ^ . 
The term in the exponent has the expansion 

E(x) = --ln, - (- + 2U--X + — *> + ... 



(7) 



% \ ( 2 4 2 6 • 7 2 



2 +2 J ^-^ 3 + --j- 
'2 4 23 2 \ /26-5 2 2 \ 

r + r + - + + - • (8) 



This result can easily be rearranged into a form 



A s = dzz- s / 2 - 2 {l-z)-^ 2 - 2 (9) 







where 

^ 1 ( 2 3 2 2 2 \ 
z{x) = -x[l --x + -x 2 + ...\ . (10) 
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For agreement with the exact result one would need to have z(l)=0.5 (the 
s and t -chanel diagrams are to cover the full range (0,1). To the present order 
in the level expansion we have z(l) ~ 0.50480 which is indeed very close to 
the exact value. A more significant observation is the fact that the main effect 
is contained in the 7 -2 factor present in the above expression . That term 
itself gives z(l) ~ 0.6. If we look up the origin of this factor, we see that it 
comes directly from the exponential higher derivative operator present in the 
vertex: since the intermediate states in the four point amplitude calculation 
are not on shell the exponential terms contribute giving the corresponding s 
dependence. 

In this nonzero momentum example, we conclude that it is advantageous 
to keep the higher derivatives exactly and that this followed by a level trunca- 
tion is likely to give a good overall approximation. Naturaly one still expects 
this expansion to only be good for certain range of momenta. 

Consider then the string field theory with the tachyon (level 0), but with 
the higher derivative terms kept exactly. The action evaluated originally in 



[pT|j reads 



C = - (d.Tf + — T 2 — Vt 3 with f = e ln ^ 2 T (11) 

The first sign that the presence of the exponential term in the cubic inter- 
action profoundly influences the nature of the problem is seen in attempting 
to evaluate the asymptotics of a possible static solution. In ordinary field 
theories, kink and lump solution can be asymptoticaly characterized as 

<t>{x) ~ 0o + aie" mx (12) 

where 0o is the constant vacuum solution and m is the physical mass of the 
scalar field. Based on this, one can write a systematic expansion for the lump 

oo 

cf>{x)= Y. a ne~ nmx (13) 

m=0 

with the classical equations providing a recursion formula for coefficients a n . 
The exponential decay e~ nmx has a physical meaning, the lump form factor 
receives a contribution from n mesons. In attempting an analogue expansion 
in the case of string field theory, one meets a surprise. After moving the 
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exponential into the kinetic term or equivalently denoting T = the string 
theory tachyon field equation reads 

((dl + l) - 2) <f>(x) = g^xf (14) 

with c = 2 In 7 = In 3 3 /4 2 . The Ansatz 

4>{x) ~ 0o + a, e"™ (15) 

leads to an eigenvalue equation 

(m 2 + l)e" c ™ 2 -2 = (16) 

for m . This equation turns out to have no real solution. With the exponen- 
tial present the nature of the lump solution has changed from that of ordinary 
field theory. In particular the decay at asymtopic infinity in string field the- 
ory has to be stronger than a simple exponential. The non-exponential decay 
of the full lump (kink) solution signals that the form factor will have a more 
complex physical meaning. This feature also necessitates a purely numerical 
approach to the problem which we attempt in the next section. 



3 Calculations and Results 

The prescence of the higher derivative terms in (14) prevent a numerical 
analysis directly in x space. Upon transforming to momentum space, one 
finds the following nonlinear integral equation 

((l _ P) e ^ - 2 ) 0(£) = g J d^ q( p(q)<P(k - q), (17) 

— * 

where k is a 25 — p dimensional vector. To solve an integral equation of this 
type, one typically discretizes the problem and solves the resulting matrix 
equation. In the case on hand, the resulting matrix equation is nonlinear and 
one is forced to search for the solution (f>(k) using a numerical minimization. 
It is computationally difficult and expensive to minimize functions with a 
large number of variables, and for this reason we have worked directly with a 
spherically symmetric ansatz for the tachyon field. Stable numerical solutions 
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were obtained using a lattice having 81 points. This implies a nonlinear 
minimization problem with 81 parameters. The choice of a good objective 
function to be used in the minimization as well as an accurate initial guess are 
crucial to obtain a nontrivial solution. Indeed, in practice, we find that the 
majority of initial guesses lead to the trivial = solution. Our objective 
function was constructed as follows: Start by making an initial guess, G(k n ) 
for the right hand side of (17). This initial guess is used to compute a value 
for the tachyon wave function 

0(U = {{l-kl)e^-2Y (18) 

The wave function (f)(k n ) can now be used to compute the value of the left 
hand side of (17). However, this involves a convolution which must be per- 
formed with care. Evaluating the convolution directly in momentum space is 
not optimal: by exploiting the spherical symmetry of the wave function, the 
convolution can reduced to the integration over a single angle and a radius. 
The integrand contains the factor (f)(p n )(f)(5 n ) with 

5 n =\k-p\ = \fp\ + k^- 2k m p n cos(0/). (19) 

In general, 5 n does not correpsond to a lattice point and one is forced to 
interpolate from the known points to find (f>(5 n ). A much more efficient way 
of performing the convolution is by transforming to x space, squaring the 
function and then transforming back to k space. By using the spherical 
symmetry of wave function, the Fourier transform can be reduced to a single 
integration. Thus, the previous double integration has been replaced by two 
single integrations, which is more efficient. In addition, the integrations only 
require a knowledge of at the lattice points. After convolving with itself 
to obtain G(k n ), the error, S to be minimized is constructed as 

£ = Y,\G(k n )-G(k n )\. (20) 

n 

The minimization was performed using the f mins .m subroutine of MATLAB, 
which employs a simplex search method. A suitable initial guess is obtained 
by taking a function which initially falls off slightly faster than a Gaussian, 
but reaches zero at some finite momentum. As an example we have shown 
the wavefunction for the D20 brane below. 
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Figl: The wave function of the D20 brane as a function of the 
radial coordinate r. The plot was obtained by taking the Fourier 
transform of the numerical solution of equation (18). 

The above form for the wave function is typical. The value of the wave 
function at the origin in position space decreased from m —0.62 for the D24 
brane to ~ —4.46 for the D18 brane. The point at which the wave function 
reaches zero is very nearly constant for all the Dp-branes considered here. In 
figure 2 we have shown the error in the D24 brane solution. The tensions of 
these solutions was evaluated directly in momentum space. For example, in 
the case of the D24 brane, we compute 

T 24 = 2vr 2 T 25 / g(l0( p )(^_ 1)eCP 2 0( _ p) + 0(p)0( _ p) 
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2n2T25 I ^(0(p)(p 2 -l)e cp2 0(-p) + 20(p)0(-p)). 
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Fig2: The error in the wave function of the D20 brane as a func- 
tion of the radial coordinate in momentum space p. The plot shows 
£/{J2n G{k n ) as a function of k n . 

In our conventions, a Dp-brane has tension 

T p = (27r) 25 - p T 25 . 
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In the table below we compare our numerically evaluated tensions, T p with 
the above values. 



V 


A T p 


Tp/Tp 


0(0) 


24 


-29.4 


0.706 


-0.63 


23 


-27.5 


0.725 


-0.86 


22 


-14.3 


0.857 


-1.19 


21 


-8.3 


0.917 


-1.63 


20 


6.4 


1.064 


-2.26 


19 


25.3 


1.253 


-3.15 


18 


64.0 


1.640 


-4.46 



Table I 



Values of the Dp brane tensions computed using the numerical 
tachyon lump solution. The parameter AT P is defined as AT P = 

Tp ^ Tp x 100. 0(0) is the value of the wavefunction at the origin in 
position space. 

Clearly there is no obstacle for the existence of lower p-Brane solutions. 
The trend is that the static lump gives systematicaly a growing tension. 
This is actually opposite to what one finds in the extreme field theory limit. 
Concerning further improvement of present results one expects (esspecialy 
for lower p) the relevance of higher massive levels. It is next important to 
study the direction of their contributions. It is also relevant to perform a 
similar study in superstring theory jl], £| f^, [16| . 
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